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The rate coefficients are calculated for trap 
loss due to excited state formation during 
5 -wave collisions of two atoms in a light 
field in a cold atomic gas near conditions 
for formation of a Bose-Einstein conden- 
sate. Blue detuning from the allowed 
atomic resonance transition causes excita- 
tion of a replusive molecular potential, 
whereas red detuning causes excitation 
when the light is tuned near a bound vibra- 
tional energy level of an attractive molecu- 
lar potential. In either case, when the light 
intensity is low and the detuning is large 
compared to the natural linewidth of the 
atomic transition, the rate coefficient for 
the coUisional loss rate is proportional to a 
molecular Franck-Condon factor. A simple 
reflection approximation formula is derived 
which permits the rate coefficient to be 
given analytically in either case. The 
Franck-Condon factor is equal to |^g(/?c)p/ 
Dc, where ^g(/?c) is the ground state 
scattering wavefunction at the Condon point 
Re, where the quasimolecule is in resonance 



with the exciting light, and Dc is the slope 
difference between ground and excited 
potentials at Rq. The analytic reflection 
approximation formula, as well as a simple 
phase- amplitude formula for the intermedi- 
ate range wavefunction, give excellent 
agreement with the results of numerical 
quantum mechanical calculations. The trap 
loss rates due to binary collisions are com- 
parable to or exceed those due to atomic 
recoil heating for a wide range of detunings 
to the blue of atomic resonance and near 
the peaks of photoassociation resonances 
for the case of red detuning. 
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1. Introduction 



Atomic collisions play a crucial role in the phenom- 
ena of Bose-Einstein condensation (BEC) [1-3]. Elastic 
momentum-transfer collisions control the rate of evapo- 
rative cooling leading to the high phase-space density 
required for the formation of a condensate. The ground 
state scattering length, also a property of elastic colli- 
sions, controls the nonlinear coupling parameter in the 
nonlinear Schrodinger equation for the condensate 
wavefunction [4] . Inelastic collisions among condensate 
atoms also determine a lower bound to the loss rate of 
the condensate (other processes may contribute to loss 
as well) and thereby help determine the condensate life- 



time. Therefore, it is crucial to understand these ground 
state collisions. Calculating their rates requires that we 
be able to compute the ground state scattering wavefunc- 
tion. Although these collisions occur in the absence of 
light, photoassociation spectroscopy offers a powerful 
probe of the ground state wavefunction [5-8]. Light also 
has the prospect for being an important probe of con- 
densate properties and for manipulating a condensate or 
even the collisions within a condensate. Therefore, it is 
crucially important to understand how light affects the 
atoms in a cold, dense atomic gas near the BEC regime. 
This article will be concerned with understanding the 
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rates of binary collisions in a weak, off-resonant light 
field in a cold atomic gas. Since only the case of a weak 
radiation field is treated, the effect of light on the binary 
collision can be described perturbatively The case of a 
strong saturating field will be treated separately [9]. The 
detuning is assumed to be large enough compared to the 
natural atomic linewidth Ja and the density A^ moderate 
enough that binary events are sufficiently isolated that 
three-body collisions can be ignored, as well as collec- 
tive effects [10]. Certainly, considering binary events is 
the first step in understanding the effect of light on 
interacting atoms. The object is to provide a simple 
framework with which to understand these binary 
events. They would need to be incorporated at a later 
stage into a more comprehensive understanding of the 
dynamics and stability of an atomic gas near BEC. 
Elsewhere, we give a brief account of the possibility of 
using photoassociation spectroscopy to probe many- 
body interactions in a condensate [11]. 

Another important reason for having a simple model 
for understanding these binary collisions in a light field 
is to understand how spectroscopic methods can be used 
as an experimental tool for probing the ground state 
scattering wavefunction. In particular, we have proposed 
that high resolution photoassociation spectroscopy 
would provide an especially sensitive probe of the 
ground state potential and scattering properties [5]. In 
fact, this technique has been used to place bounds on 
scattering lengths for Li [8], Rb [6], and Na [12, 13] 
ground state collisions. Even tighter bounds should be 
possible with improved experiments. In particular, we 
will show below how the spectroscopy of very cold 
gases might be use to measure properties of the ground 
state wavefunction. 

The primary effect of light on cold trapped atoms is 
to cause loss processes by giving the atoms enough 
energy to escape the trap. For example, a single atom 
scatters light of frequency o) detuned Aa = oj-oja from 
the atomic resonance frequency oja at the rate. 



/atomic — '^X A 



(1) 



where /2a = ^/IttIIc d^ is the Rabi frequency at laser 
intensity / for the atomic transition with transition 
dipole 6?A, and we assume \Ap\ » /2a and \Ap\ » Ja, Aa 
is defined so as to be positive for blue detuning and 
negative for red. Since each scattering event gives a 
recoil momentum kick of hk^ = hwlc to the atom, and 
since the thermal momentum of the atom is comparable 
to, or smaller than, hki^, this atomic light scattering 
heats the atoms, causing their loss from the trap. 



Binary events also cause loss of atoms from a trap, 
but with a rate which is very different for the cases of 
red and blue detuning. Figure 1 illustrates the two cases 
schematically, using the molecular potential energy 
curves for the quasimolecule comprised of the two col- 
liding atoms. Blue detuning excites the ground state 
atoms to a repulsive molecular state, causing the produc- 
tion of excited atoms with kinetic energy Ee = Eg-\- ^Aa. 
This heating has been discussed in the context of optical 
shielding [14,15,9]. The probability for this free-free 
process is a smooth function of detuning. On the other 
hand, red detuning excites the ground state atoms to a 
bound eigenstate of the attractive molecular state. This 
only happens with high probability if the light is tuned 
close to exact resonance with the energy E^ of the bound 
vibrational state v. Since the decay of the bound state 
mostly leads to products that are not trapped, this pro- 
cess gives rise to a trap-loss photoassociation spectrum, 
which is now well studied for the alkali dimers [16-21]. 
The probability for this free-bound process is a sharply 
peaked function of detuning. We have every reason to 
expect that a gas near BEC conditions, or even a conden- 
sate itself, will possess a detailed photoassociation spec- 
trum for detuning on the red side of Wa- 

When the laser intensity / is sufficiently low, the 
probability for both red and blue detuning loss processes 
can be described by a Franck-Condon factor between 
the ground and the excited state radial wavefunctions. 
For both cases, we will give here a very simple analytic 
expression for the Franck-Condon factor that is in excel- 
lent agreement with fully quantal calculations. We will 
show that the Franck-Condon factor is proportional to 
I ^g(/?c, E)\'^, the square of the ground state wavefunction 
at the Condon point Re, where the quasimolecule is in 
exact resonance with the light. At Re the difference 
between excited and ground potentials exactly matches 
^(o . Our expression for the binary collision rate exhibits 
the proper quantum threshold law behavior as £" ^ 0. It 
predicts that the binary collisional loss rate per atom, 
Tbinary, also varlcs as A~a and at typical BEC densities 
generally is comparable to or exceeds the atomic scatter- 
ing rate, Eq. (1), for a wide range of blue detuning and 
greatly exceeds the atomic scatering rate when (o is 
tuned near a photoassociation resonance. 

We will first describe the model and the basic proper- 
ties of the ground state wavefunction. Then we will 
discuss the binary collision rates for the case of blue 
detuning and the reflection approximation for the free- 
free Franck-Condon factors. Then we will describe the 
modifications needed to describe the case of red 
detuning and free-bound Franck-Condon factors. 
Finally, we will consider some of the implications for a 
cold atomic gas near BEC. 
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Quasimolecular Absorption 



Potential 

Energy 




Intemuclear separation R 

Fig. 1. Schematic diagram of quasimolecular potential energy 
curves for two atoms with an allowed S ^ P transition, showing 
free-bound and free-free absorption for respective red and blue detun- 
ing from the atomic resonance frequency oja- 

2. Model 

The binary collision rates for the two cases described 
in Fig. 1 will be expressed in terms of Franck-Condon 
factors between the ground and excited state wavefunc- 
tions. However, it is very desirable to formulate the 
problem using the field-dressed molecular picture of a 
collision in a light field. This procedure is well- 
documented elsewhere [22-27], even for cold collisions 
[28-30,5,14], and need not be described in detail here. 
A two-state model with a single ground state and a 
single molecular excited state will be assumed. In the 
weak-field limit it is straightforward to generalize to 
multichannel ground and excited states for real alkali 
systems with hyperfine structure. In fact, we have now 
developed computer codes that let us include both 
ground and excited state molecular hyperfine structure 
for alkali species for optical transitions in the weak field 
limit [12, 13]. The ground state potential is Vg(/?) and 
the excited state potential is V^iR )- ^ Aa, where both Vg 
and Ve ^ as /? ^ 00. The ground and excited states are 
coupled optically by the interaction matrix element 
^fl{R), where the molecular Rabi frequency is 



dressed potentials cross at the Condon point Rq, the 
point of quasimolecular resonance, where in the un- 
dressed picture of Fig. 1 , the difference of upper and 
lower potential curves matches the photon energy /tcj: 



hiOA + Ve(/?c) - Vg(/?c) = fiio 



V.{Rc)-VJRc) = Ma. 



(3) 



(4) 



In this paper we will carry out purely model illustra- 
tive calculations. For this purpose it is sufficient to use 
analytic Lennard-Jones potentials for the ground and 
excited states, 

i = g or e and n = 6 for the ground state van der Waals 
potential and /i = 3 for the excited state resonant dipole 
potential. The long range variation is C^/R^ and C^/R^ 
respectively. We use a reduced mass characteristic of 
sodium. It is difficult to show on one graph the poten- 
tials over the very large range of energy and distance 
involved in cold collisions. Figure 2 shows on a logarith- 
mic scale the magnitude |y(/?)|/^B (where ^b is the 
Boltzmann constant) of the long range potentials for 
the ground and excited states of the Na2 molecule, for 
which we take Ce = 1500 atomic units (e^al, where 
e = electron charge and Gq = Bohr radius = 0. 0529 nm) 
and IC3I = 10 atomic units (e^al). Figure 3 shows on a 
linear scale the ground and attractive and repulsive 
excited state potentials, as well as illustrative wavefunc- 
tions for each. Since the excited potential is orders of 
magnitude stronger than the ground state at long range, 
it is an excellent approximation to neglect the ground 
state contribution in calculating the Condon point: 






=^ Rc = 



(6) 



Figure 2 also indicates the detunings associated with 
each R , taken as a Condon point. Note that the ground 
state interaction strength is not negligible at intermediate 
/?, but is stronger than 1 |ulK whenever R is less than 
about 300 Gq. 



n{R) = V277//C d{R) = b{R)nA 



(2) 



where d(R) is the molecular transition dipole, and 
^ b{R) ^ 2/V3 relates the molecular and atomic 
Rabi frequencies. The particular value of b depends on 
the molecular states involved. The ground and excited 



3. Ground State Wavefunction 
3.1 Long Range Form 

We need the ground state wavefunction to use in the 
reflection formula for the Franck-Condon factors 
we derive below. When the collision energy E is low 
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Fig. 2. Magnitude of long range ground state van der Waals and 
excited state resonant dipole interaction potentials, shown using log 
scales. The interaction strengths were taken to be Ce = 1500 atomic 
units and C3 = 10 atomic units respectively, characteristic of Na atom 
interactions. The scale on the right hand side of the figure gives the 
detuning A^ corresponding to the Condon points associated with the 
resonant dipole curve. 
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Fig. 3. Long range ground and attractive and repulsive potentials on 
a linear scale for the same long range parameters as for Fig. 2. Typical 
free and bound excited state wavefunctions for a Na reduced mass are 
also shown, for an exit channel with 5 cm"' kinetic energy and for a 
bound state with a binding energy of 4.87 cm" ' . The 2 fjiK ground state 
wavefunction is nearly linear between 60 a^ and 200 a^. 



enough, only collisions with zero relative angular mo- 
mentum of nuclear motion contribute to scattering cross 
sections. The ^-wave ground state wavefunction takes on 
the following well-known behavior as /? — > 00: 



\'%{E)) - d'^B 



~H^T 



sm(k„R + Tjg) 



(7) 



where the Vz/I/tt^^ factor ensures energy normal- 
ization, |<^(£)|^(£'))p = 8(E-E'). The asymptotic 
kinetic energy is £' = /t^kLUfji^ where jjl = 171^12 is the 
reduced mass, the deBroglie wavelength is Aa= = lir/ko,, 
and the velocity v^ = ^ klfx. It is sometimes more 
convenient to write the energy normalization factor as 






(huo)" 



(8) 



If E is small enough, the asymptotic phase r]^ = nir- 
^gAs, where As is the scattering length and n is the 
number of bound states supported by the ground state 
potential; the nir term can be ignored, since it only 
introduces an inessential phase factor. The scattering 
length is a property of the whole potential and for real 
atoms is very sensitive to small uncertainties in the short 
range, chemical bonding, range ofR. However, the scat- 
tering length provides a good parameter to characterize 
the effect on the long range wavefunction of the com- 
plete potential. 

Figure 4a illustrates typical ground state wavefunc- 
tions (in their real form without the trivial complex 
phase factor e^"^-) with positive, zero, and negative scat- 
tering lengths. These were calculated by making slight 
changes in the model potential parameters. The inner 
region wavefunction, say for R <30 «o, is characterized 
by rapid oscillations as the local kinetic energy 
increases due to the acceleration by the attractive ground 
state potential. The intermediate range wavefunction, 
say for SO ao < R < XJA, is nearly linear in R and 
extrapolates to an intercept near R = A^: 



^JR, E) - e 



■- &)' 



/2 k^{R-Ad 



(9) 



A much more accurate representation of the ground 
state wavefunction in the intermediate range is found by 
correcting for the variation in phase and amplitude 
due to the nonzero ground state potential. The details 
of this derivation, based on an approximate treatment of 
the rigorous Milne equation for the phase-amplitude 
form of the wavefunction, are given in Appendix A. The 
result is that a more correct form of the long range 
wavefunction is found from Eqs. (60), (69), and (72) 



^(/?, E) 






X a(R) sm(k^p(R)X 



where 



.<.)=. -(I)- 



(10) 



(11) 
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Fig. 4. (a) Comparison of characteristic ground state near-threshold 
wavefunctions at 1.4 |jlK collision energy. The wavefunctions were 
calculated using a Na reduced mass and three different hypothetical 
model potentials which yield scattering lengths A^ = - 96 ao, 0, uq, and 
+ 101 ao. The intermediate range wavefunction projects to a node near 
R = As. Of course, the scattering length for any real atomic species is 
fixed by the actual interaction potential and can not be varied (except 
possibly by introducing external fields), (b) Comparison of the scatter- 
ing wavefunction for the potential for which A^ = -I- 101 ao with the last 
bound state wavefunction (bound state binding energy is 2.6 mK), 
when both are given a common WKB normalization near the distance 
of the deepest part of the ground state potential. 



MR>'R-A.-i (iy«. 



(12) 



Equation (10) only applies in the region R » 7?b so that 
the correction term (Rb/R T « 1 , where 



and 120 ao for a model Na potential with As = + 84.3 Uq. 
Such a magnitude of As is near the actual value for 
doubly spin-polarized Na [31,32,12]. The actual wave- 
function node at 87.5 ao is shifted 3.2 ao further out than 
As, as accurately predicted by Eqs. (10) and (74). A 
much larger shift of 1 8 ao in node position is predicted 
for ''Rb. 
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Fig. 5. (a) Comparison of the exact calculated ground state wave- 
function (points) with the predictions of the phase amplitude approx- 
imation, Eq. (10), using a Na reduced mass and a potential for which 
As = 84.3 uq. The actual node position atR = 87. 5 ao is shifted to larger 
R . The node position is independent of collision energy, as predicted 
byEq. (12). 



3.2 Short Range Form 

Before leaving this discussion of the ground state 
wavefunction, it is useful to comment on the nature of 
the short range wavefunction as £^ — > 0. When the colli- 
sion energy E decreases towards zero, there will always 
be a characteristic quantum threshold connection be- 
tween the asymptotic and short range wavefunction. It is 
possible to define a characteristic distance Rq where the 
WKB connection between the asymptotic and short 
range wavefunctions strongly breaks down for collision 
energies E < Eq [33, 34]. A necessary condition for 
onset of the threshold law behavior is E « Eq. The 
values of Rq and Eq are determined by requiring that the 
maximum in the function d\(R, E)/dR satisfy the 
following condition: 



^B = 



10^^ 



(13) 



dXiRQ^Eg) ^ 1 

dR 2' 



(14) 



is a constant independent of energy with dimensionality 
of distance. It has a value of 42 Uq for Na and 77 Uq for 
Rb. Figure 5 compares the exact wavefunction and the 
result of Eq. (10) in the near linear region between 60 Uq 



where the local deBroglie wavelength X(R,E) is defined 
by Eq. (63) in the Appendix A. For all E < Eq there is 
a range of distances around7?2wheredA(7?,£^)/d7? > 1/2 
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and the WKB approximation fails. Using Eqs. (17) and 
(18) of [34], 



Rn = 



20' 



2"'(7a)' 



7?B= 1. 0437?B 



(15) 



1 /i4\3 n\\!i 
where a = — \ — \ ( - 1 = 7.0419. Thus, we find that 

Rb = 0.96Rq, independent of potential parameters (since 
both have identical scalings with Ce and fx), and the 
defined parameters Rb and Rq turn out to be fortuitously 
very close in magnitude. Since each of these only de- 
fines a qualitative range of distance associated with a 
change in the wavefunction, they can be used inter- 
changably. The WKB breakdown region is actually very 
broad (An example for a He mass is in Fig. 3 of 
Ref. [33]), and becomes broader as energy is lowered 
below Eq= 19 mK for Na or ^^ = 1 ■ 5 mK for Rb [34]. 
Therefore, significant departure from the WKB form 
begins at distances well to the right of Rq, and Eqs. (10) 
and (11) show that the amplitude of the intermediate 
range wavefunction remains much closer to its asymp- 
totic amplit ude g<= l/ vk^ than to the local WKB ampli- 
tude cc l/\/k(R, E). Thus, the amplitude of the ground 
state wavefunction does not appreciably exhibit the ac- 
celeration due to the ground state potential unless R is on 
the order of or less than Rq ~ Rb. This is one of the 
reasons why a semiclassical analysis fails in the long 
range region as £" ^ 0. 

The form of the wavefunction in Eq. (60) is com- 
pletely rigorous at all distances and energies, providing 
that the Milne equation described in Appendix A is 
solved exactly. References [33, 34] show that the ground 
state wavefunction for the inner region R « Rq is well 
approximated by the form 



^l(R,E) = (A,k^y'^c' 



JllL 



\R.E) 



X sin()S7^^7?,£)) 



= (A,k^y'^ ^J''^^R,E=Ol 



(16) 



(17) 



function near threshold and the wavefuncton for the last 
bound state, both normalized at short range to the same 
WKB form. Scattering wavefunctions near threshold, as 
well as the last bound state wavefunction, when given 
this common normalization, are almost indistinguish- 
able for R < Rq. The physical reason for this is that the 
local amplitude and phase are determined by the ground 
state potential in a region where it is very deep com- 
pared to the initial collision energy or bound state bind- 
ing energy. 

Both the inner range and outer range approximate 
wavefunctions, Eqs. (16) and (10), are proportional to 
vk^, in accordance with the threshold law requirement 
that the ground state wavefunction for R « \lka> every- 
where vanishes as v /tco. This ensures that the probabil- 
ity for inelastic processes due to the ground state en- 
trance channel vanish as /:« as £" — > 0. This is the basic 
threshold law for inelastic collision rates which must be 
satisfied for light-induced inelastic events for either red 
or blue detuning. 



4. Collisions for Blue Detuning 

The case of cold collisions in a blue detuned light field 
has been extensively studied recently, both experimen- 
tally and theoretically [14,15,9,35^1]. Figure 4a shows 
the ground and excited state potentials for the case of a 
very large detuning of 5 cm"\ or 300 GHz, from atomic 
resonance. If the light intensity is strong enough (this 
will be defined more precisely below), the ground state 
atoms approaching one another encounter an effective 
repulsive interaction near the Condon point Rq and are 
repelled. This gives rise to the phenomenon of optical 
shielding [15,9,35^1], by which the light prevents the 
atoms from approaching closer than Rq, thereby reduc- 
ing the rate of ground state processes which require the 
atoms to be closer together than Rq. The light also mod- 
ifies the ground state elastic scattering rate and scatter- 
ing length, and results in the creation of hot excited state 
atoms which share the kinetic energy released, hAp, 
[14,15,9]. Since we are considering the weak field limit 
here, the shielding effect is small and we will concen- 
trate on the latter energy release process. The event rate 
coefficient for this process is 



where A, is an energy-independent constant having t he 
dimension of length and aJ^(R, E) = l/Vk^ (R, E) is 
the WKB amplitude. Equation (16) shows that 
the overall shape of the inner range wavefunction, 
i^wKB + ^^ ^^ ^^^j^ -g independent of energy, whereas the 
overall amplitude decreases as {Aik^y^ as £" — > 0. This 
is illustrated in Fig. 4b, which shows a scattering wave- 



K^{A^) = (^ P^{E, A^)j = ^^ P^iE, A^)}, (18) 



where the brackets imply an average over the distribu- 
tion of ground state velocities Ug, and 



P^(E, A^) = |5eg(£, A^f 



(19) 
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is the probability for the event by which two ground 
state atoms collide in a light field and produce one 
excited and one ground state atom, both of which have 
enough kinetic energy to escape any weak trapping 
potential. 

When the radiative coupling is small, the radiative 
distorted wave approximation [24,26,27], can be used 
for the 5 -matrix element: 

(20) 

where Veg(^) = ^0(R) from Eq. (2). The kinetic energy 
in the asymptotic ground and excited state channels are 



E and E + ^Aa respectively The asymptotic excited state 
wavefunction is: 

\W:(E, Aa)\) - ^~'^'{-^;^^y' sin(/:e,./? + 7/e) . 

(21) 

The ground and excited state wavefunctions are also 
illustrated in Fig. 6a. The excited state wavefunction 
oscillates much more rapidly than the ground state one 
because of the much greater kinetic energy in the exit 
channel. 

By assuming that the radiative coupling Veg(^) is 
either independent of Ror slowly varying with R near 
Re, then this term can be removed from the integrand in 
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Fig. 6. (a) Ground and repulsive excited state field-dressed potentials for a blue detuning Ap, of 
5 cm"^ = 300 GHz and laser intensity 7=0. The crossing at /? = /?c becomes an avoided crossing 
when 1^0. Ground and excited state wavefunctions are also shown, (b) Similar figure for the case 
of a red detuning Aa of -4.87 cm~\ in resonance with the v = 84 level of the attractive excited 
potential. 
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Eq. (20), giving the radiative distorted wave result that 
Pe(£, ^a) = Att^VI K^(£, ^a)|^(£))|^ (22) 
= A7T^VlF,,{E,A^\ (23) 

where Vc = V^^iRc) and 

Feg(£, A^) = K^(E, A^)\^;(E))\' (24) 

is a free-free Franck-Condon factor. Although direct 
evaluation of F^^ by numerical quadrature does not 
converge, since both wavefunctions oscillate with finite 
amplitude to /? = oo, there are several methods available 
for evaluating F^^ numerically [27], as well as approxi- 
mate formulas based on the local nature of the integrand 
near the Condon point Rq. Here we use a two-state close 
coupled scattering calculation of the collision in a weak 
radiation field with an /? -independent Rabi frequency, 
and extract S^^ie) from the asymptotic field-dressed 
wavefunction [24-27]. Then the Franck-Condon factor, 
independent of the laser intensity chosen for the calcula- 
tion, is calculated from: 



F.UE, A^) 



K(E, A^ 



(25) 



F!,^(E, A^) 



hvcdc 



where 



Dr 



(V.(R)-VJR)) 



dR 



_3C3_ J_ I . I 



(26) 

(27) 
(28) 



is the difference in potential slop es evaluated at Re and 
V(. = Vg(Rc) = \/2{E - V^{Rq)Iil is the local velocity at 
Re. The second equation follows from neglecting V^ in 
comparison to K- The corresponding semiclassical 
treatment in the field-dressed picture of Fig. 6a is given 
by the Landau-Zener (LZ) approximation [14,15,9] for 



Pt\E, A^) = 2 e-^^ (1 - e-^^) - 2 Ac- 



(29) 



The latter approximation applies when the dimension- 
less LZ adiabaticity parameter 

(30) 



hvcDc 



where S^^{E, A^) is calculated using a standard close 
coupling scattering code. This method can be general- 
ized to a multichannel version when there are more than 
two channels [27], and thus could be used for atoms 
with hyperfine structure [12, 13]. 

Figure 7 shows F^^ calculated for collision energy 
Elk^ = 2 /aK using a model ground state potential with 
the reduced mass for Na collisions and a calculated 
scattering length of +84.3 a^. The detuning range corre- 
sponding to /? = 20 (3o to 200 (3o is 8 THz to 8 GHz. The 
Franck-Condon pattern shows an oscillatory pattern 
versus Re or detuning. 



5. Reflection Approximation 

We will now turn our attention to developing approx- 
imations that will give considerable physical insight into 
the nature of light scattering at very low collision energy. 
Let us first look at semiclassical approximations. We 
will assume a single Condon point Re. Thus, our analy- 
sis does not apply to states with significant contributions 
from two Condon points, such as the low vibrational 
levels of the Og state of Na2 [12, 13]. A semiclassical 
approximation for the free-free Franck-Condon factor 
F^^(E) is given by Eq. (46) of the reference [42] (we 
remove the uninteresting Vc factor): 



is small compared to unity. This condition also ensures 
the low intensity regime of linear variation of P^ with 
light intenstiy. Using Eq. (25) to obtain Fi^(E) from 
Eq. (29) gives exactly the same expression as in 
Eq. (26). Note that both Eqs. (26) and (30) have left out 
the phase factor sm\l3.(Rc) - PdRc))- Actually, the 
phase factor has been replaced by its average value, 1/2. 
Note that P^^ given by Eq. (29), even with the phase 
factor, does not satisfy the threshold law requirement to 
vary as ^a= as £" ^ 0. In fact, F^^(E) from Eq. (26) is in 
very poor agreement with the quantum mechanical 
/eg(£') in Fig. 7 and is not shown in the figure. 

The usual stationary phase derivation of Eq. (26) can 
be adapted to give the correct threshold law as £" ^ 0. 
The details are given in Appendix B. The result is a key 
result of this paper, a simple reflection formula: 



F^,(E,A^) = ^J^;(Rc,Ef. 



(31) 



The Franck-Condon factor is directly proportional to the 
square of the ground state wavefunction. It only depends 
on the excited state through the slope term Dc. This 
expression explicitly satisfies the threshold law, since at 
all distances R « Aoc, | ^gp is proportional to ^a= when E 
is sufficiently smaller than Eq [33, 34]. When R > Rq, 
the quantum wavefunction is well-approximated by 
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Fig. 7. Comparison of exact and approximate evaluations of the Franck-Condon factors 
for different Condon points Rq. The correspondance of detuning to Re is indicated on the 
top horizontal axis. The free-free Franck-Condon factors have dimensionality of 
(energy)"^, where here energy is expressed in cm"^ units. Since the same magnitude of the 
C3 coefficient was use for the repulsive and attractive curves in our model calculations, the 
exact free-bound Franck-Condon factors are seen to lie on nearly the same curve as the 
free- free ones, when the former are divided by dEJdv to normalize them per unit energy. 
This is as predicted by Eqs. (31) and (41). The reflection approximation formula of Eq. 
(31), |^g(/?c)|^/£*c, is in excellent agreement with the exact calculations over the whole 
range shown. The intermediate range phase- amplitude formula, Eq. (32), is good agree- 
ment for /? > 70 ao. (a) logarithmic scale (b) linear scale showing region near outer nodes. 
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Eqs. (10) -(12), and the intermediate range phase- 
amplitude approximation to F^^ is: 



the whole range of temperatures encountered for 
trapped laser cooled atoms, including the E — > limit. 



FltiE. ^a) : 



hv^.Dc 



alsm\k:^pc). 



(32) 



6. Collisions for Red Detuning 



where ac = a(Rc) and pc = p(Rc)- The basic difference 
between this formula and the Landau-Zener one, 
Eq. (26), other than the phase factor, is the presence of 
the asymptotic ground state velocity in the denominator 
instead of the local velocity tt- If ^c is in the range 
Rq « Re « Aoo/2 so that the sin function is linear in its 
argument, then 






(33) 



The equivalent approximation for P^ is found from 
Eq. (23), 



P.iE, AjO 






(34) 



We can carry out a very similar treatment for the case 
of the free-bound Franck-Condon factor for red detun- 
ing. Figure 6b illustrates the dressed potentials and 
wavefunctions for a red detuning case. The basic differ- 
ence from blue detuning is that the attractive potential 
supports discrete bound states. When the laser fre- 
quency (o is tuned near the position of a particular bound 
level with vibrational quantum number v, the level can 
be excited and decay to some product p , which is no 
longer trapped. The decay process is primarily sponta- 
neous emission leading to formation of molecular spe- 
cies or hot atoms [28, 42]. In the field dressed picture of 
Fig. 6b, the bound state is a scattering resonance embed- 
ded in the continuum very close to threshold. The rate 
coefficient for red detuning is also given by Eq. ( 1 8), but 
the 5 -matrix element is given by a modified resonant 
scattering form described by Napolitano et al. [5]: 



Figure 7 also shows the predictions of the F^ and F^^ 
formulas compared to the exact results for 2 |ulK colli- 
sions. The full reflection formula is an excellent approx- 
imation over a wide range of Condon points, and even 
remains very good into about 20 Uq where the ground 
state has accelerated the atoms to a velocity much higher 
than their initial one. The F^^ approximation formula is 
also excellent at intermediate range, until R begins to 
approach Rb- Rq. Figure 8 compares the exact and F^^ 
results for other collision energies up to 2 mK. We see 
that the reflection approximation is very good for almost 



ypjsiE, V, Aa) 



(35) 



^ 



where A^ is the detuning relative to the position By of the 
bound state. The total decay rate of the excited bound 
state vis % = yp -H jsiE, v, Aa) + %, where jp is the rate 
by which the bound state resonance decays to the de- 
tected product, jsiE, V, Aa) is the stimulated emission 
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Fig. 8. Comparison of the exact free-free Franck Condon factors and the predictions of 
the phase- amplitude formula, Eq. (32), for higher collision energies, E/kB = 2 |jlK, 200 |ulK, 
and 2000 fjiK = 2 fjiK. The breakdown of the reflection approximation is seen at the larger 
Re for 2000 |ulK. 
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rate back to the ground state continuum, and y^ is the 
decay rate due to any other undetected processes (such 
as molecular predissociation). The critical element for 
our purposes is the ys(E, v, Aa) factor, which describes 
how the cold colliding atoms couple to the excited bound 
state. For low light intensity, this factor is given by the 
Fermi golden rule expression [5]: 



y,(E, u, Aa) 



277 I 



[^Xv)\V.,(R)\^,(E))\\ (36) 



where ^eCu) is the unit-normalized bound state wave- 
function of the excited bound state. Using the same 
factorization approximation as for Eq. (23), this can be 
written in terms of a bound-free Franck-Condon factor. 



277 



%(£, V, /1a) = 



vh\{%(v)\n(E)w 



^ Vh F,,(E, V, A^). 



(37) 



(38) 



The assumption is that Vc is evaluated at a single 
Condon point; if there is more than one Condon point, 
the interference between the different amplitudes from 
each Condon point would need to be taken into account. 
The derivation of the reflection approximation 
formula for the bound- free Franck-Condon factor F^^{E, 
V, Z\a) follows immediately from the same treatment in 
Appendix B that we used in deriving the free-free 
Franck-Condon factor. The essential difference is that 
the unit-normalized bound state wavefunction can also 
be written in energy- normalized phase-amplitude form 
just like Eq. (16), as is commonly done in generalized 
multichannel quantum defect theory [44]. It is only 
necessary to introduce the vibrational spacing 



dE^ 
dv 



= hv^. 



(39) 



where Vy is the vibrational frequency, the number of 
complete vibrational cycles per unit time. Then the unit 
normalized bound state wavefunction has the form [44] 



^.{R.v) = 



dE^y^ 
dv, 



2fji 



a,(R,v)sm(l3.(R,v)y 

(40) 



where a and ^ can be determined quantum mechani- 
cally from Eqs. (61) and (62), or alternatively, from the 
WKB form, Eq. (64). The derivation in Appendix A, 
adapted for a right hand turning point, carries through to 
give 



Feg(£, u, Aa) =^^^ \n(Rc. E)\\ (41) 

where ^(^c, E) can be calculated exactly or taken 
from Eq. (10)) above. The vibrational spacing function 
can be evaluated approximately from the discrete level 
spacing, 

dEy _ Ey+l - Ey-l 

dv^ 2 



(42) 



or more rigorously from the wavefunction at some short 
range point ^o: 



dv ' 



2(1 



{ko\^. 



(Ro, vf 



]_ 

ko 



d^,(Ro, v) 



dR 



')■ 



(43) 



where ^o = K(Ro)- This equation immediately follows 
from Eq. (40) upon using the WKB forms in Eqs. 
(77)-(78). The result depends only weakly on the choice 
of Ro as long as it is in the classical region of the 
potential well. We use Eq. (43) in the numerical results 
described below. 

Although the free-free and free-bound Franck- 
Condon factors have different dimensionality, they may 
be compared by evaluating the free-bound Franck- 
Condon factor per unit energy, or F^^(v, AA)(dEy /dv)~\ 
Figure 7 also shows the results of evaluating this func- 
tion for a number of discrete bound levels using an 
attractive potential that has the same magnitude of the 
C3 long range potential coefficient as for the repulsive 
state (see Fig. 3). Each bound level defines a distinct 
Condon point nearly equal to the outer classical turning 
point. The free-bound and free-free Franck-Condon 
factors fall on almost exactly the same curve in Fig. 7, 
even when Re is as small as 20 Uq. This is expected from 
the fact that the free-free and free-bound transitions in 
the model calculation have the same ground state wave- 
function and Dc at the same Condon points. The results 
in Fig. 7 show that the reflection approximation is just as 
good for bound states as for free states. 

Now we can easily derive from Eqs. (38) and (41) a 
simple expression for the stimulated decay rate of the 
excited bound state to the ground continuum, ys(E, Aa): 



r.(£,^A)=^^ hVy\^,(Rc, E)\' 

27T Vl , 4 2-2/7 ^ 



={ 



\6Tr^Vl 2,, 2,,, 



(44) 

(45) 

(46) 
(47) 
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The probability Pe{E) in Eq. (47), defined by the ex- 
pression in braces in Eq. (46), is identical in form to the 
probability in Eq. (34) for a free-free transition. The 
expression in Eq. (47) thus has a pleasing physical inter- 
pretation: the total rate of decay out of the bound state 
is the probability of decay during one vibration cycle 
(corresponding to a "complete" collision with incoming 
and outgoing parts) times the frequency of vibration 
(number of cycles per second). The overall probability 
for the photoassociation transition, given by the reso- 
nance scattering expression in Eq. (35), explicitly 
exhibits the threshold law form due to the js factor in the 
numerator. 

The form in Eq. (44), proportional to the square of the 
ground state wavefunction, will also give the threshold 
law form for/?- and 6f-waves, etc. It also explicitly shows 
why the photoassociation spectra map out the nodal 
structure of the ground state wavefunction in the inten- 
sity pattern versus vibrational quantum number [28, 17]: 
as detuning Ap, changes from level to level, the 
Franck-Condon factor tracks the ground state wavefunc- 
tion at the changing /?c- This effect is expected to be 
much cleaner in very low temperature spectra, where 
only ^ -waves contribute, than in most of the existing 
higher temperature experiments, where s-,p-, d-wavcs, 
etc., all may contribute. The p- and d-wavcs do not 
generally have the same nodal structure in the interme- 
diate range as the ^-wave does. 



This expression is independent of collision temperature, 
as is should be. Using Eqs. (2), (6) and (28) and the fact 
that we can write [43, 42] 



= /3%a(^)'. 



(50) 



the rate coefficient can be algebraically transformed to 
the form 



2K^'^%A^) 



2b^ 

377 



Ta 



m ^/=- 



(51) 



In Eq. (50) /s is a fraction of order unity depending on 
the particular molecular states involved (In Hund's case 
(a) it is 1/2 for a 2 state and 1/4 for a 11 state; in general 
it will be a function of R that can be worked out for any 
particular state). The following definitions are used in 
Eq. (51): 



be = b(Rc) 



Aa 



(52) 



(53) 



/c 



^'-mi'-k-im- <«) 



7. Collision Rates Near BEC 

In order to estimate the effect of binary collisions in 
a cold dense atomic gas, it is necessary to compare the 
rate for the binary process with the atomic scattering 
rate yatomic in Eq. (1). The binary event rate is 



Tbinary = ^e(^A)A^, 



(48) 



Recall that b(R) is defined by Eq. (2). The Ib^^BTT 
factor in Eq. (51) is a dimensionless molecular physics 
factor depending on the specific transitions involved, 
and it has a maximum value of 

977 r 



where K^ is given by Eq. (18). For the purpose of esti- 
mating the relative importance of atomic and binary 
light scattering events in causing loss of trapped atoms 
due to heating, it is necessary to compare yatomic and 
27binary, slncc two atoms are lost per binary event. For 
the large detuning and moderate density conditions we 
consider, the influence of any collective effects on light 
scattering is not likely to alter significantly these basic 
magnitudes. We will base our estimates on the interme- 
diate range formulas for R > Rb, Eqs. (34) and (46), 
since these apply to a wide range of detuning from near 
resonance to beyond 100 GHz. 

Using Eqs. (34) in (18), the rate coefficient for trap 
loss events for blue detuning in the intermediate range is 



2^^"^(ziA) = 



hDc 



2 2 

acPc- 



(49) 



The density A^a corresponds to one atom per cubic Aa. 
The density near BEC conditions in a cold atomic gas is 
typically larger than A^a ■ The factor /c has a nominal 
order of magnitude unity and reflects the phase of the 
ground state wavefunction. It will have a node near As if 
As is positive and in the intermediate range. 
Using Eq. (51), we find 



2vu-"^ 1 A^ 

-^ /binary J_ ^v 

/atomic ' ^^X 



(55) 



This fraction is nearly independent of the blue detuning 
Aa over a wide range, since/c ~ 1 (except near the node). 
Since Nx = 5xm'^ cm "' for Na and 2X10" cm"' for Rb, 
we see that the order of magnitude of the light-induced 
binary collision loss rate is comparable to or larger than 
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the atomic light scattering rate when the density is in the 
rangeof 10" to 10^^ cm'^ 

It is also possible to get a simple expression for the 
rate coefficient for trap loss due to a bound state photo- 
association resonance by using Eqs. (47) and (35) in 
Eq. (18). The thermal average can be greatly simplified 
at very low T when k^T is very small compared to the 
linewidth y^ of the bound state photoassociation reso- 
nance V. Since the numerator of the resonance scatter- 
ing 5 -matrix element, Eq. (35), satisfies the threshold 
law form, the energy-dependence disappears in the 
numerator of the resonance scatttering rate coefficient 
expression. Then we are left only with the energy depen- 
dence in the E-Ay term in the denominator of Eq. (35). 
But since k^T is small compared to either of the terms 
in the denominator, the E dependence is negligible and 
can be dropped. We also assume y^ = y^, that is, each 
decay results in a loss event, since the stimulated decay 
rate is assumed to be small {y^ « y^). Thus, the thermal 
average rate coefficient due to resonance v becomes, in 
the r ^ limit, 



minimum rate coefficient from summimg over all such 
contributions is ~ 2Ki^V(4A^^). Thus, we estimate 



/ tTUo 



^«» = lir'"<"'''''j(iFTw5) 



Jp 



= 2K^J- 



^vjp 



(Ayf + (yJ2f 



(56) 



(57) 



In Eq. (57) 



2Ki^\A^)=^^^alpl 



(58) 



is identical in form to 2^e^"^ from Eq. (49). It is the 
single-cycle loss rate coefficient, due to passing the 
crossing once in each direction. The net rate is deter- 
mined by the resonance expression. We can envision two 
possible limits. The first is the exact resonance case 
when the laser is tuned to the bound level so Ay = 0. In 
this case, the resonance enhancement factor over the 
single pass result reduces to 4vyTy, where r^ = l/y^, is the 
decay time of level v. But v^Tv is just the number A^^ of 
complete cycles in one lifetime. So the maximum loss 
rate coefficient is 2^^(i;) = 8^e ^ A^... The second limit 
is when the detuning is as far as possible from reso- 
nance. Assuming equal spacing of vibrational levels, the 
largest detuning is just half of the vibrational spacing; 
thus, the maximum detuning is just A^ = ttVv » y^ This 
gives \l{7T^Ny) « 1 for the factor multiplying IK^^ in 
Eq. (57). The minimum value the rate coefficient can 
take is 2^^" = 2K^^^I(7t^N^). Of course, we must add the 
off-resonance contributions of all levels v -\, v + 1 , 
etc., to get the total of /-resonant rate coefficient. If the 
vibrational spacing is assumed to be uniform, the 



1 



{216:^) < {2ICf) < (4N.)(2Ki'^). 



(59) 



From Eq. (59) and (55), and the fact that A^^ » 1 
whenever Z\a » y^, we see that 2yb1nary will be a sharply 
peaked function that is much smaller than 2ybin^y (for 
the same \Aa\) when o) is tuned between resonances, but 
which is much larger than 2ybin^y when the o) is tuned 
to a photoassociation resonance. Therefore, by detuning 
to the red side of atomic resonance, the effect of binary 
collisions can be made to be either very large or very 
small compared to yatomic- The cold atomic gas will have 
a prominent photoassociation spectrum, and the rate for 
photodestruction at the peak positions will measure the 
ground state pair correlation function oc|^g(7?c)P at the 
Condon points for the transitions. There is a discrete set 
of Condon points Rc(y), corresponding to the outer turn- 
ing points for the vibrational levels of the upper molec- 
ular state. There are several well-known alkali dimer 
transitions that have already been studied [7] in much 
hotter traps than the < 1 |uiK traps that are now available. 
The intriguing question remains whether 3 -body or col- 
lective effects might perturb the binary photoassociation 
spectrum, causing broadening or shifts in the resonance 
line profiles. Kagan et al. [45] have suggested the possi- 
bility of observing quantum correlations in optical prop- 
erties of low temperature spin-polarized gases. The 
photoassociation spectrum of a cold dense atomic gas 
near BEC might be a very fruitful source of experimen- 
tal and theoretical investigation of 2- and 3-body interac- 
tions in such systems [11]. 

Certainly the analysis in this paper shows that light 
can induce significant loss processes, in comparison to 
that caused by free atom light scattering, due to colli- 
sions between atoms in a cold atomic gas near 1 |uiK or 
below at densities > 10'^ cm"^. If the density should 
increase to > 10'^ cm~^, the collisional loss rate will 
dominate that due to atomic recoil heating, except 
between photoassociation resonances for large enough 
red detuning. For the case of the MIT Na trap [3] the 
blue detuning of the laser used to plug the hole in the 
magnetic trap was so far off-resonance that any atomic 
or binary light scattering did not prevent the formation 
of a condensate. But for much smaller detunings the 
simple Franck-Condon formulas we have derived here 
give a way to estimate the binary collision rates for trap 
loss via excited state production. More detailed treat- 
ment, including incorporating the actual quasimolecular 
structure of alkali dimer states, could be readily incor- 
porated into this framework. We have also developed a 
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theory for treating the nonperturbative effect of hght, 
including multiphoton effects and the modification of 
the ground state scattering length [9]. Careful attention 
to collective effects also needs to be given for the case 
of higher densities. In any case, it seems very unlikely 
that a condensate could be stable for very long in the 
presence of light tuned within a few hundred GHz of 
atomic resonance. 



8. Appendix A. Milne Equation Analysis 
of Asymptotic Wavefunction 

The use of the exact Milne quantum mechanical equa- 
tions for the phase and amplitude of the wavefunction 
lets us do an exact analysis of the asymptotic wavefunc- 
tion. Let the wavefunction be writen in phase-amplitude 
energy- normalized form: 



W(R, E) = [^j c^(R, E) sin(^(/?, E)X (60) 



where a satisfies the Milne equation [46, 47], 



dR 



-^ k\R,E)] a(R,E) 



1 



a(R,Ey 



0. (61) 



The phase j8 is calculated from a by 



and the local wave vector is 



(62) 



^(^'^^ = ^=Vl^(^-^(^» 



(63) 



The asymptotic potential is V(R) = - Ce/R^ 

In the semiclassical picture, we have the identifica- 
tion of a with the WKB amplitude 



a^^(/?,£) = 



1 



\/k(R,E) 



(64) 



and the above forms are equivalent to the WKB wave- 
function. The exact solution to the quantum equations 
[Eqs. (61) and (62)] give an exact quantum wavefunc- 
tion. Here we will use an asymptotic analysis for the 
wavefunction at moderately long range. 

At long range the amplitude becomes constant, and 
very large as £^ ^ 0: 



1 



V£ 



2^' 



(65) 



Therefore, we want to calculate the departure of this 
large, nearly independent of R, amplitude from its 
asymptotic limit. We do this by approximating the last 
term of Eq. (61) by 1/a^ ~ kl a. Then Eq. (61) becomes 



^+(k\R)-kJ)a = 



dR 



j-^V(R)a = 



dR^ 



{^C,)j,a. 



(66) 
(67) 

(68) 



Equation (68) can be immediately integrated to give 






(69) 

(70) 



where R^ is defined by Eq. (13). When R = 20''^/?b = 
2.1 Rb, the amplitude a{R) has only changed by 5 % 
from its asymptotic value aa>. When R = Rb, the correc- 
tion term is unity, and the asymptotic approximation 
used in setting up the approximate Eq. (66) breaks 
down. Therefore, the amplitude expression in Eq. (70) 
can only be applied for R » Rb- 

Now we can turn our attention to deriving the correc- 
tion to the asymptotic phase. This is done by substituting 
the above amplitude, Eq. (70), into the phase [Eq. (62)]. 
The asymptotic phase at some suitable large value of 
R = Rc^ (a distance well beyond the range of the poten- 
tial) is j8a> = ko.(R:>o - As), where A^ is the 5 -wave scatter- 
inglength. The phase function at R is calculated by inte- 
grating backwards from /?«=, 



l3(R) = k4Ro.-As)-^ 

= k4Ro.-As)-k 






t{Rf 
dR 



(71) 



» k^iR^-A,) -k^j (l + 2 [^y)<iR 



= UR-A, 



(72) 
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Combining the corrected phase and amplitude 
expressions, Eqs. (72) and (70), in Eq. (60), we obtain 
the expression given in Eq. (10) in the text. Note that the 
normalization stays close to its asymptotic value, in spite 
of the fact that the local k(R,E) can be much larger than 
^00 in that part of the moderately long range region 
where |yg| » E (see Fig. 2). Thus, the semiclassical 
local WKB normalization does not apply. Note also that 
the node of the wavefunction at R^ near As for the posi- 
tive scattering length case is shifted from its value 
projected from the asymptotic A^. This is because the 
phase near A^ has not yet reached its asymptotic limiting 
value. This shift in the actual node to larger R can 
readily be estimated from Eq. (72): 



2fji 



Rn-As 



2 //?i 



3 \R> 



and Rn ~ As implies /?„ ~ As[ 1 + :r( -r- 



Rn = 



2 Ri 



(73) 



(74) 



where 



feg(£, 4a) = ^ |/eg(£, 4a)P, 

/eg(£, Aa) = I agOeSinySgSin^Sed/? 
Jo 

~ - I ftgOe COS()3g - )3e)d/?. 



(79) 



(80) 



(81) 



In the last line we have made the standard approximation 
of neglecting the rapidly oscillating phase sum term. 
The phase difference is expanded in a Taylor series 
about a point R^: 

cj>(R, E) = p,(R, E) - p,(R, £e) (82) 

= Z70 + b,(R-R,) + i b2(R-R,f + . . . , (83) 



where 



9. Appendix B. Reflection Formula 

This derivation of the reflection approximation ex- 
pression follows Jablonski [48] and the Appendix of Ref. 
[49]. We begin by writing the ground and excited state 
wavefunctions in energy-normalized phase-amplitude 
form 

%(R. E) = i^Y a,(R, E) sin(^g(/?, E)) (75) 



^e(/? , £e) = l^] "^ a,(R , E, ) sm(l3.(R , E, )). (76) 



where E^ = E -\- ^Aa. The ground state amplitude and 
phase are given by Eqs. (70) and (72), whereas the 
excited state quantities can be taken to be the standard 
WKB ones: 



a,(R,E,) = 



Vk(r7e7) 



I3XR,E,)= f ^e(/?',£e)d/?' + ^77. 



(77) 



(78) 



where Rt is the inner turning point of the repulsive 
potential. The Franck-Condon factor is 



bo = )3g(/?p. 


E)-p,(R„E,) 


(84) 


= )3g(/?p,£)-| 


(85) 


'" dR 


dp^iR, £.) 
Rp dR 


Rp 


1 


Y - URp, E) 


(87) 


" a{Rp, E 


'-k.-K{Rf,E,) 


(88) 


u d%{R,E) 
"'- dR' 


d%{R, £.) 


(89) 


"'' dR^ 



dK(R,E,) 
dR 



fJiPc 

^X(Rp) 



Rp 



(90) 



(91) 



In the normal semiclassical theory, Rp is taken to be the 
Condon point, where k^ = k^ so that bi = 0. Here we 
simply take the point where bi vanishes. Since the 
excited state potential is so steep with respect to the 
effectively flat ground state potential, it is an excellent 
approximation to take /?p = /?t = ^c to all be equivalent. 
Then from Eq. (78), P^Rp, E) = 77/4. The expression for 
b2 follows from neglecting the ground state contribution 
in comparison to the excited state one. 
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We thus evaluate 



\L 



UE) = aJRc, E)a,iR, E,) - cosibo + ^ x')dx 






= ^2^] a,(Rc, E) sm(l3,(Rc, E)). (92) 



£ 



where we have used | cos( bo ± -^ x^]dx 



^cos(z,o±f). 



(93) 



Using Eqs. (75) and (92) in Eq. (79) gives the desired 
reflection approximation to the Franck-Condon factor: 






(94) 
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